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We study the possible decay of a coherently oscillating scalar field, interpreted as dark matter, 
into light fermions. Specifically, we consider a scalar field with sub-eV mass decaying into a Fermi 
sea of neutrinos. We recognize the similarity between our scenario and inflationary preheating where 
a coherently oscillating scalar field decays into standard model particles. Like the case of fermionic 
preheating, we find that Pauli blocking controls the dark matter decay into the neutrino sea. The 
radius of the Fermi sphere depends on the expansion of the universe leading to a time varying 
equation of state of dark matter. This makes the scenario very rich and we show that the decay rate 
might be different at different cosmological epochs. We categorize this in two interesting regimes 
and then study the cosmological perturbations to find the impact on structure formation. We find 
that the decay may help alleviating some of the standard problems related to cold dark matter. 

PACS numbers: 95.35.+d, 98.80.-k 



I. INTRODUCTION 

Dark matter has become an extremely interesting area 
of research in both cosmology and particle physics. From 
the particle physics point of view it can be thermal 
WIMPs, axions, Kaluza Klein states, etc. Though su- 
persymmetric (SUSY) models predict its mass to be of 
the order the electro-weak scale, there are viable models 
of dark matter where its mass can be as low as sub-eV, 
for example axion-like dark matter. Especially the direct 
and indirect search for dark matter has narrowed down 
the parameter space of these well studied candidates to 
a large extent. Hence it is highly likely that dark mat- 
ter is of much more exotic nature than thought of. In 
addition, there may be a requirement for more compli- 
cated physics such as interactions with dark energy or 
with neutrinos. Similarity between the neutrino mass 
and the present dark energy density scale has inspired 
people to look for a connection between the two. Now 
we know that the normal active neutrino cannot be a vi- 
able dark matter candidate because of its free-streaming 
ability. But the existence of sub-eV neutrino mass might 
point towards richer sub-eV scale physics. In fact there 
has been a few interesting works p], 0] where new states 
of sub-eV masses are present in the dark sector. The 
reason is, if the dark sector interacts only gravitation- 
ally with the standard model sector, a TeV scale SUSY 
breaking in SM would predict scalar particles of mass 
T , e X ~ 10 _3 eV in the dark sector. Scalar dark matter 
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of milli eV mass with a possible coupling to neutrinos has 
been discussed in Q . Also moduli of sub-eV mass can 
easily arise from string compactification 

If in nature dark matter arises from such a low energy 
scale, we would expect it to decay into light fermions 
like neutrinos through a Yukawa type of coupling. GeV 
scale dark matter decay (and annihilation) to neutri- 
nos has drawn lots of recent interest [oT-[ll|. especially 



in the context of recent cosmic ra y m easurements and 
the DAM A/LIBRA experiment (l2Ttl3j. But in our case, 
dark matter is of sub-eV mass and the signature of its 
decay into neutrinos is mainly cosmological, especially in 
structure formation. Recently in a different context, the 
non-thermal wimp miracle [14j was introduced where a 
scalar of TeV mass decays into a stable dark matter par- 
ticle. So the decay of a scalar particle can lead to very 
rich phenomenology in cosmological context. 

The possible decay of scalar dark matter into light neu- 
trinos is an effect which could potentially help to under- 
stand the apparent surplus of power on small scales in 
simulations containing normal CDM This sur- 

plus could for instance be an indication that CDM is 
simply clustering too much on small scales and we need 
some mechanism to reduce their gravitational interac- 
tion. This is where the decay could play a role - see [l8[ 
for a similar idea. In this paper we study the nature of 
the decay and its possible signature in structure forma- 
tion. As we consider an axion-like scalar dark matter, 
it is undergoing a coherent oscillation and might decay 
into neutrinos through a parametric excitation. In that 
case, the process will have many similarities with infla- 
tionary fermionic preheating |19l423 |. Therefore we dub 
the process Preheating Dark Matter. 

Cosmology with decaying/interactin g, da rk matter has 
been an interesting topic of research [23^29] in recent 
times as it gives a probe to detect dark matter indi- 
rectly through its effects on structure formation. On 
top of that, if it decays into dark energy it might give 
rise to a unified description of dark matter and dark 
energy. In most of the studies [3fJ, l3l| , the rate of en- 
ergy transfer from dark matter to other components (eg. 
dark energy, radiation or neutrinos) has been empirically 
assumed driven by mathematical simplicity. Here we 
present a concrete model of dark matter decay to light 
fermionic states like neutrinos and then study its imprint 
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on structure formation. In particular, we derive the de- 
cay rate of dark matter as a function of redshift using 
the theory of fermionic preheating. We find that in the 
initial stage the decay rate is faster and determined by 
parametric resonance. But at late times the parametric 
production ceases and redshifts of fermionic modes con- 
trol the decay rate. This time varying decay rate makes 
the phenomenology rich and offers a prominent imprint 
on the structure formation, something which might be 
experimentally probed by near future experiments. 

The plan of the paper is as follows: In section II we 
discuss the particle physics aspect of our scenario. In 
section III, we incorporate the idea of inflationary pre- 
heating for scalar dark matter decaying into neutrinos. 
In section IV, we identify different epochs of decay and 
derive the background evolution, i.e. how dark matter 
and neutrino energy densities evolve in presence of the 
decay. In section V, we derive the perturbation equa- 
tions for our scenario and in section VI we obtain tem- 
perature anisotropy and matter power spectra. Finally, 
we summarize our work in section VII. 



II. PHENOMENOLOGICAL MODEL 

Here we discuss how our scenario fits into a particle 
physics set up. Scalar fields of sub-eV masses are com- 
mon in different particle physics models. In many models 
of TeV scale gauge mediated supersymmetry breaking, a 
gravitationally coupled dark sector contains a scalar of 
sub-eV mass. Also, string compactifications generically 
predict axion like scalars such as the dilaton and large 
numbers of moduli [32J, |33[ whose mass can easily be in 
the sub-eV range. 

As we are interested in a sub-eV scalar field which cou- 
ples to light fermions like the standard model neutrino, 
our set up is inspired by the models of mass varying neu- 
trino dark energy [34j where a sub-eV mass scalar couples 
to the standard model neutrino - see also [35l - l39j . Espe- 
cially, we refer to a model of supersymmctric neutrino 
dark energy, where multiple scalars of such low mass are 
present and can couple to neutrinos. In such a theory, it 
has been shown that the scalar potential takes the form 
of the well known hybrid inflation potential. We refer 
readers to |l| for details and will briefly discuss here. 
The Lagrangian for mass varying neutrino dark energy is 
given by 

C D m D vN + kANN + h.c. + V(A) (1) 

where mo is the Dirac mass and v is the left chiral Weyl 
field representing active neutrinos. N is the righthanded 
heavy fermion and A is the scalar, k is some Yukawa 
coupling. 

In supersymmetric models of neutrino dark energy 
where v,N,A is promoted to supcrficlds l,n,a, the su- 
perpotential takes the form 

W = nann + rn^ln (2) 



After taking quantum corrections into account, it has 
been shown that this leads to a hybrid inflation kind of 
potential. Depending on the temperature of the universe, 
the scalar either remains trapped at a metastable mini- 
mum playing the role of dark energy or it rolls off and 
starts oscillating coherently, behaving like cold dark mat- 
ter. Following a simple model @, the Lagrangian looks 
like 

C =An 2 ?/>3 + 2An 2 V'2^3 + "13^3^3 + mi^)%vi 

+ V S usy + V sof t +V e + h.C. (3) 

where 

V susy = 4X 2 n 2 n 2 3 + \ 2 n\, (4) 

and 

V soft = m\n\ - m\n\ + 03713 (5) 

The terms in V e are included in order to generate a Ma- 
jorana mass for the neutrino in the vacuum. In this kind 
of theory the superpartner sneutrinos (here denoted by 
^3) can easily be of sub-eV mass and play the role of 
dark matter. Also it can easily couple to light fermions 
(neutrinos). From now on, we will switch our focus to 
cosmological effects of such a model. 

III. PREHEATING FROM SCALAR DARK 
MATTER 

We are essentially interested in a light scalar field dark 
matter of sub eV mass which has a coupling to an ultra 
light fermion which for our case is the standard model 
neutrino. We consider decays of such dark matter into 
neutrinos, though our framework is true for decay into 
any fermion. Following the mechanism of inflationary 
preheating, in this section we will understand the physi- 
cal nature of the decay and will clarify different regimes 
of decay. We will see that the decay rate changes as the 
universe expands due to time evolution of a resonance 
parameter which controls the parametric excitation of 
neutrinos. 

We mainly follow the work [22j on fermionic preheat- 
ing and apply that to our scenario. So we refer to this 
work for detailed derivation of the equations. Briefly, to 
find the number density of created fermions through the 
preheating mechanism in an expanding background, one 
derives a mode equation using the original Dirac equation 
with the Friedmann-Lemaitre-Robertson- Walker metric. 
It has been shown that the comoving number density of 
created fermions can be obtained by solving for a mode 
function Xk(t). For a Yukawa type coupling X^tpip, the 
mode equation is given by 

x" k + [k 2 + (m + ^f) 2 - iyfqf']X k = 0, (6) 

where 4>of{t) is the background solution for the time 
evolution of the oscillating scalar field, k = is the 
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my, 



and the reso- 



dimensionless fermion mass, m 
a 2 <a 2 

nance parameter q = — ^ . These three parameters com- 

<p 

pletely determine the parametric production of fermions. 
We consider the oscillation of the field with the usual 
quadratic potential V = ^m 2 (j> 2 as this is a good ap- 
proximation around minima of any potential. The term 
(m + y/qf) can be thought of as an effective mass of the 
fermion. As the scalar field oscillates, the effective mass 
itself will oscillate around zero and the parametric pro- 
duction of fermions is enhanced when the effective mass 
crosses zero. It has been shown numerically that nk(t) 
oscillates and due to Pauli blocking its maximum value 
never crosses unity. But for decay into bosonic particles 
it is not bounded by unity. We stress that the behavior of 
this parametric production is considerably different than 
the perturbative decay process <f> — > ipip where the decay 
rate is given by T ~ . 

In the above mode equation, expansion of the universe 
has been neglected which may only be true at very late 
times where the Hubble parameter drops. To get a full 
understanding, one must include the expansion of the 
universe. This alters two aspects. The parameters q 
and k now become time dependent. More specifically we 



get q 



A-^(t) 



and the physical momentum p 



where a(t) is the scale factor of the universe. As a re- 
sult, the periodic modulation of the comoving number 
density does not hold anymore. For large values of the 
resonance parameter (q > 1), the calculation of para- 
metric production becomes, in fact, simple. Luckily, we 
will see later that for our case q ^> 1 for large periods 
of (cosmological) time. Using the method of successive 
scattering for fermions, it has been shown that due 
to the loss of periodicity of nk the production of fermions 
happens through a stochastic filling of a Fermi sphere up 
to a Fermi radius kf which depends on scale factor a(t) 
and is given by 



(7) 



Now to find the exact dependence, one needs to know 
how q(t) changes with scale factor. As we are using a 
quadratic potential for the scalar, the solution for the 
scalar field for this case is well known. Oscillation of <fi 
in this case is given by the asymptotic solution 



a 3/2 

Using this it is easy to derive 



cos(t). 



1/4 1/4 



(8) 



(9) 



where qo 

So, as the universe expands, the Fermi sphere also ex- 
pands producing more and more neutrinos. But the res- 
onance parameter decreases as the amplitude of oscilla- 
tion drops due to Hubble friction and at some point the 



Fermi sphere stops expanding when q(t) becomes of the 
order of unity. At this regime, the redshift of fermionic 
modes due to Hubble expansion is fast enough to prevent 
the parametric excitation. Finally, fermions will be pro- 
duced with a much lower rate in the perturbative regime 
and perturbative processes continue unless < 2m,/,. 



IV. DIFFERENT DECAY REGIMES 

Using the above results now we can focus on produc- 
tion of neutrinos and its time evolution. Here we assume 
that parametrically produced neutrinos mix with other 
relic neutrinos and acquire the same temperature through 
thermalization. It is instructive to note that as the mass 
of the scalar is way less than in the usual inflationary 
preheating scenario, we would get parametric excitation 
until very late times. From the previous discussion, we 
have learnt that the resonance parameter q is very cru- 
cial to determine the nature of the decay and q itself is 
time dependent. Now we will discuss the two different 
decay regimes and the transition time between them for 
our simple model with a quadratic potential. 



A. Regime I: Expanding Fermi radius (q >> 1) 

During early stages of parametric production 

q{t )^^l^2X 2 ^. (10) 



As we are interested in scalar mass of the order of ~ 
1CP 3 eV, almost all over the cosmic history until today, 
Pdm > (10 _3 eU) 4 . Now if the coupling constant A is 
of the order of unity, we still get parametric production 
at very late times. But for smaller couplings parametric 
excitation stops at earlier redshift when q(t) ~ 1. Later 
we will take different choices of the coupling A and study 
how it affects the formation of structure. The produced 
neutrino number density is obtained through the volume 
of the Fermi sphere with radius 



= q-> -a'' - x a 1 = g 1/4 a - 3/4 , (11) 
This can be used to calculate the neu- 



wherc q = — 

* rn 

trino density 



Pv 



F d 3 fc = 8vrA 2 im^ 2 = 8ttA 2 pdm, (12) 



where we have used pdm = ^ rn ^,4' ■ It is important to 
note that the neutrino energy density is proportional to 
the local dark matter density. Using this fact and the 
continuity equation for the total dark matter and neu- 
trino fluid, we can find the evolution of the dark matter 
energy density and hence neutrino energy density. The 
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continuity equation for the dark matter and neutrino as 
a whole reads 



ptot + 3-ffptot (1 + tfltot) = 0, 
I77' Pto* = Pv + Pbm and P to 



(13) 



P, J + P^ 



DM- 



where w tot = 
Eq. [T3] can be split up into the two components 

Pdm + 3Hp DM = -Q (14) 

and 



Pv 



4Hp u = Q, 



(15) 



where Q represents the decay rate from dark matter to 
neutrinos and we have taken advantage of the fact that 
Pdm = and P v = \p v . Combining Eqs.[TJ[Tl and [TJ 
we get the relations 



Pdm = Pdm 



(?)' 



p v = 8ttX pdm, 



(16) 



where the i denotes the value at some fixed time(e.g. 

today) and i = 3 1 t f 3 8 2 ^ ■ So here we clearly see that 
due to parametric production, dark matter no longer red- 
shifts as 1/a 3 . Its effective equation of state changes from 
zero to slightly higher values. The higher the coupling 
A, the higher the deviation. This effective equation of 
state which corresponds to the value we would get if we 
did not know about the coupling between dark matter 
and neutrinos can be calculated from a revised version of 
Eq. EH PDM + 3Ppdm(1 + w cS ) = 0. The result is 



Mefi 



1. 



(17) 



We note that for A — > 0, it gives the right limit for the 
equation of state u) e g — > 0. 



B. Regime II: Fermi radius stops expanding (q ~ 1) 

In the second regime, parametric excitations weaken 
due to the drop in resonance parameter q. During this 
regime, q ~ 1 and the radius of the physical Fermi sphere 
has approached the constant value kp ~ m^. This means 
the decay is controlled by the redshifts of Fermi momen- 
tum due to the expansion. As the universe expands, the 
Fermi momentum drops, opening up space in the Fermi 
sphere. This space is immediately filled up by the scalar 
field decaying into neutrinos. The regime may be impor- 
tant for structure formation if the decay into neutrinos 
can cause a substantial decrease in the dark matter den- 
sity. This is possible when pdm ~ because, in this 
case, decay of each DM particle to a neutrino causes a 
significant decrease in the dark matter energy density. 
As the dark matter mass is of the order sub eV in our 
model, this can happen only at late times. This can en- 
hance the late ISW effect thus modifying the structure 



formation on large scales. Again, using Eq. [12] we easily 
obtain p„ = 47rm^. From the continuity equations Eq. 1141 
and Eq. [T5] we get the relations 



Pdm = 

7 

Pu = J, 



(18) 



where 7 = 167rm^. So we see that the neutrino energy 
density is constant in this regime, only the dark matter 
density dilutes. This is what we expect, because if there 
were no decay into neutrinos, its density would simply 
redshift like the standard model neutrino governed by 
the Hubble expansion. But here, as soon as phase space 
opens in the Fermi sea of neutrinos due to cooling of the 
universe, it gets refilled by the decay from dark matter 
thus keeping its density constant. 

In this regime the effective equation of state for dark 
matter can be calculated to be 



Wcff 



Pdm 



(19) 



which depending on model parameters can deviate an 
appreciable amount toward the present - see Fig. [T] 
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FIG. 1: The effective equation of state of dark matter in the 
preheating dark matter scenario with the following parameter 
choices: Dotted green: JIdm.o = 0.21, A = 0.001 with Q^o = 
0.022, dot-dashed blue: ft D M,o = 0.06, A = 0.05 with Sl v , = 
0.081, and dashed red: SIdm.o = 0.12, A = 0.03 with fl v ,o = 
0.068. 



Transition redshift 



In this subsection we find the transition redshift be- 
tween the two epochs in terms of model parameters 
and A. Patching the two regimes together at a scale fac- 
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tor ax we arrive at the relations 
p v = 8ttA 2 pdm 

PDM = 2P 



Pv = \ 
PDM = ' 



{pdm 



5) (^) £ 



for a < ax 



for a > ax 



(20) 

where the denotes present day values, and ax can be 
determined from the preheating dark matter parameters 
as 



ax 



(pdm + j) 

2A 7 ~ 3 



1/3 



(21) 



V. PERTURBATION ANALYSIS 

In order to study the implications of preheating dark 
matter, we perform a cosmological perturbation analysis 
in the synchronous gauge in which a line element is given 
by 



ds 2 = a 2 (r) [-cZt 2 + (f/y + h ij )dx i dx j ] 



(22) 



where i,j = 1,2, 3, ?yy is the Minkowski space metric, Ziy- 
is the perturbation to the metric and we are using comov- 
ing coordinates x^ = (x, r) in a spatially flat background 
space-time. We follow the procedure given by Ref. [40l j . 
in which, to linear order in the perturbations, the stress- 
energy tensor is given by 



n 



-(p + Sp) 



T? = (p + P) Vl = ~T* 

Ti = {P + SP)S) + Y^, £* = 0, 



(23) 



where the perturbations to energy density and pressure 
are defined as Sp = p — p and SP = P — P, is the 
anisotropic shear perturbation, and Vi is the coordinate 
velocity of the fluid 1 . The latter is a small quantity and 
can be treated as a perturbation of the same order as 
dp and SP. Instead of working with the velocity itself 
we use the divergence defined as 9 — ik l Vi. Similarly, 
instead of the anisotropic shear perturbation, we use the 

shear stress a. This is defined as a = Still — iL^iiEi. 

(p+p) 

The conservation of energy and momentum for our 
coupled fluid implies that the covariant derivative of the 
stress-energy tensor is 0. 



1 \H fluid — "f 1 



1 a/3 1 + 1 a/3 1 



0. 



(24) 



1 For further information about cosmological perturbation theory 
see Ref.[47l 



However, for the individual components in the fluid it is 
slightly different 



-pfj,a 

1 ;ti CDM 



-SQ 
SQ, 



(25) 



where the 8Q can be determined directly from Eq. [TJ] 
Using the time-time(OO) component of the stress-energy 
tensor from Eq. [23] and combining with Eq. [25] above, 
we get the equation of motion for the individual density 
contrasts Si = Using the space-space(ii) components 
in the same way will give us the time evolution of 9i . 

Hence we arrive at the equations of motion for the DM 
and neutrino components for q 1 



S v 



'CDM 




#CDM — 0. 

Similarly, for q ~ 1 we get 



(26) 




'CDM 



£ PCDM 



1 



% = -AH9 U + k 2 -5 V - a v 



'CDM 



0. 



(27) 



It is the effect of the term containing A and the term 
containing the neutrino energy density p v in the evolution 
of the CDM density contrast that separates the evolution 
of perturbations in the preheating dark matter case from 
the normal case, where the decay of CDM is not permit- 
ted. 



VI. RESULTS 

In this section we present the numerical results we 
obtained by using the equations from the two previ- 
ous sections. We modified the publicly available CMB- 
FAST code [4l| to include preheating dark matter to get 
both the matter power spectrum and the temperature 
anisotropy spectrum. This code is developed to calculate 
the linear CMB anisotropy spectra based on integration 
over the sources along the photon past light cone, but also 
outputs transfer functions from which the linear matter 
power spectrum can be calculated. 

In our analysis we keep the epoch of matter-radiation 
equality fixed, and the only free parameters are the cur- 
rent value of !]dm,o and the parameter A. In addition, 
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we keep the amount of baryons today fixed at Qb = 0.05 
and choose the normalized Hubble expansion rate at the 
value /iHubbic = 0.7. We include one species of massless 
neutrinos produced in the decay as well as the three stan- 
dard model neutrinos, which for simplicity are assumed 
to have a degenerate mass spectrum with m v = 1.5x 10~ 3 
eV. We assume that the neutrinos produced in the decay 
mix with standard neutrinos, although this assumption 
makes no qualitative difference to the results. 

The impact of preheating dark matter on the temper- 
ature anisotropy spectrum, can be seen in Fig. [2] The 
most apparent difference from the spectrum of ACDM 
can be seen on largest angular scales, I < 100 (correspond- 
ing roughly to a degree), although for some choices of pa- 
rameters the positions and relative heights of the peaks 
are also affected. We generally observe an increase in 
power on scales, 20 < I < 100, whereas on scales Z < 10 
we see an increase or a decrease in power depending 
on the model parameters. For scales Z < 100, the dom- 
inant contribution to the CMB temperature anisotropy 
spectrum comes from the Integrated Sachs- Wolfe effect, 
which arises because of the evolution of the gravitational 
potentials encountered by the photon on its journey from 
the last-scattering surface. The modification to the cos- 
mological background because of preheating dark matter 
can be quite significant, as we saw in the last section. In 
particular in the q ~ 1 regime, where the neutrino energy 
density becomes constant, which leads to a second term 
in the evolution of the CDM perturbation in Eq. [27J 

On the largest angular scales I < 20, the dominant con- 
tribution to the CMB temperature anisotropy spectrum 
comes from the late-time Integrated Sachs- Wolfe effect 
(IS W) . This effect is a result of the universe entering an 
epoch of rapid expansion as it becomes dominated by 
dark energy. In this epoch, photons moving into gravi- 
tational potential wells will get a boost as the potential 
well is decaying and becomes slightly shallower while the 
photon is passing through it - and vice versa for gravita- 
tional hills. It is clear that this effect depends intimately 
on many different parameters such as Odm.o and n u ,o- In 
our model, dark matter is being transformed into light 
neutrinos - i.e. radiation - most efficiently at late times, 
and hence we expect an effect on the largest scales. In the 
context of interacting dark matter-dark energy models, 
the ISW effect has been studied in 

Turning our attention to the matter power spectra, we 
know that the linear matter power spectrum is extremely 
well determined by SDSS (see Ref. ^) and WMAP (see 
Ref. [ifl) on intermediate and large scales. And in ad- 
dition, Ly-a forest data have some constraints on the 
small scales - see e.g. (4(|. Hence we normalize our mat- 
ter power spectra such that they coincide with matter 
power spectra obtained from using normal ACDM at the 
largest scales (these are also the latest to have entered 
the horizon). The results are presented in Fig. [3] for dif- 
ferent values of A and ilDM,o- A small damping on small 
angular scales seems to be generic, similar to standard 
models of CDM and hot dark matter, where a similar 
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FIG. 2: The temperature anisotropy spectrum as a function 
of the angular modes in the preheating dark matter scenario 
with the following parameter choices: Solid black: ACDM 
with £7dm,o = 0.24, J1a,o = 0.71, dotted green: JIdm.o = 0.21, 
A = 0.001 with 0^,0 = 0.022, dot-dashed blue: f2 D M,o = 0.06, 
A = 0.05 with Q„ t o = 0.081, and dashed red: floM, o = 0.12, 
A = 0.03 with Q„ o = 0.068. 



reduction in power is achieved. We note that in order to 
comply with e.g. supernova data (Ref. (42j ) we cannot 
change Om.o drastically. Of course we do turn CDM into 
neutrinos that redshift as radiation - only slightly faster 
than CDM. Hence we have some room to change ^m,o 
and still be in agreement with data. 

The matter power spectra for the different parameter 
choices are agreeing relatively well with ACDM as re- 
sult of the small A-value. Still, we do notice the small 
reduction on the smallest scales which can be probed by 
CMBFAST. This is to be expected since part of the CDM 
responsible for the gravitational wells is decaying into 
neutrinos which undergo free-streaming on these small 
scales. This prevents them from clustering and since we 
are creating an appreciable amount of neutrinos due to 
the decay of CDM, we generally expect this reduction of 
power on small scales. 

On even smaller scales we also expect a considerable 
effect, which will reduce the clustering ability of CDM. 
Unfortunately, we cannot probe those scales satisfacto- 
rily with CMBFAST - as it is using linear perturbation 
analysis and we expect the evolution to be highly non- 
linear. The investigation of the effect of preheating dark 
matter on non-linear scales is beyond the scope of the 
present work and will be postponed to a future paper. 



VII. CONCLUSION 

In this work, we have considered a coherently oscillat- 
ing scalar field of sub-eV mass which behaves as dark 
matter. Due to its coupling, it slowly decays into neu- 
trinos as the universe expands until today. The decay 
rate as a function of redshift has been derived follow- 
ing the physics of inflationary preheating of a scalar into 
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FIG. 3: Top: The matter power spectra as a function 
of k[Mpc _1 ] in the preheating dark matter scenario with 
the following parameter choices: Solid black: ACDM with 
Hdm.o = 0.24, Oa.o = 0.71, dotted green: floM, o = 0.21, 
A = 0.001 with fi„,o = 0.022, dot-dashed blue: f2 D M,o = 0.06, 
A = 0.05 with ft^o = 0.081, and dashed red: S7dm,o = 
0.12, A = 0.03 with il u ,o = 0.068. Bottom: The differ- 
ence between the preheating dark matter model and ACDM 
(5 2 (ACDM)/<5 2 (PHDM)) as a function of k. 



fermions. We find that the decay rate is modulated by 
Pauli blocking and the expansion of the universe, giving 
us rich physics of dark matter decay into a neutrino sea. 

We studied the effect of the decay on structure for- 
mation and obtained spectra for the anisotropics in the 



cosmic microwave background and matter power spectra. 
For the parameters proposed in this paper, we showed 
that given the decay we are able to slightly reduce the 
amount of power on small scales in the matter power 
spectra - something which seems to be required from data 
- while at the same time being in good agreement with 
SDSS and WMAP observations. 

Interestingly, the proposed scenario leads to features in 
the temperature anisotropy spectrum - which can be seen 
as a prominent late ISW effect and slight modifications 
to the second and third peaks. Consequently, as future 
direction, it would be interesting to study the late ISW 
effect in detail, from which we will be able to constrain 
the model parameters more effectively. In addition, we 
would like to do a follow-up GOSMOMG analysis using 
the newest data available as well as to examine the effect 
of preheating dark matter on structure formation on non- 
linear scales. After such an analysis it will be more clear 
what the best choice of preheating dark matter parame- 
ters is such that a comparison with future Planck data, 
for instance, will be easier. Furthermore we expect fu- 
ture weak lcnsing surveys to be useful in constraining our 
scenario since they can provide insight into the statistics 
of the dark matter distribution - hence they can (hope- 
fully) shed some light on what happens on small (and 
large) scales of gravitational clustering. 
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